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ABSTRACT 

We discuss the possibilities for primordial black holes (PBHs) to grow via the accretion of dark 
matter. In agreement with previous works, we find that accretion during the radiation-dominated era 
does not lead to a significant mass increase. However, during matter-domination, PBHs may grow by 
up to two orders of magnitude in mass through the acquisition of large dark matter halos. We discuss 
the possibility of PBHs being an important component in dark matter halos of galaxies as well as 
their potential to explain the ultra-luminous x-ray sources (ULXs) observed in nearby galactic disks. 
We point out that although PBHs are ruled out as the dominant component of dark matter, there 
is still a great deal of parameter space open to them playing a role in the modern-day universe. For 
example, a primordial halo population of PBHs each at 10 2 ' 5 M Q making up 0.1% of the dark matter 
grow to 10 4 ' 5 M© via the accumulation of dark matter halos to account for ~ 10% of the dark matter 
mass by a redshift of z w 30. These intermediate mass black holes may then "light up" when passing 
through molecular clouds, becoming visible as ULXs at the present day, or they may form the seeds 
for supermassive black holes at the centers of galaxies. 

Subject headings: black hole physics — accretion — dark matter — galaxies: formation 
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1. INTRODUCTION 

It is well-established that supermassive black holes 
(SMBHs) with masses in the range itibh ~ 10 6 — 
10 9 - 5 M^ reside in t he centers of spheroidal systems 
l|Bernardi et al.ll200lh One can make a convincing case 
that these have grown largely through accretion, with the 
consequent energy emission observed in ele ctromagnetic 
output and jets at an efficiency of e ~ 0.1 (SoltaiJ ri982t 
lYu fc Tremainel 12002(1 . Observations of distant quasars 
have shown us that these SMBHs are already in place 
by redshifts of 6 and greater, but the mechanism of their 
formation remains a mystery. 

Motivated by this question, we examine the potential 
for primordial black holes (PBHs) to grow through ac- 
cretion to become seed masses for SMBHs. Primordial 
black holes, defined as black holes forming in the early 
un iverse without stellar progen itors, were first propo sed 
bv lZel'dovich fc Novikovl ljl967j) and IHawkinel l)l971|) as 
a possible consequence of the extremely high densities 
achieved in the Big Bang model. If they do in deed form 
in the early universe and can avoid evaporation ijHawkind 
11975(1 up to the present day, they must still exist and they 
may be important. We use a combination of analytical 
and numerical methods to follow PBH growth through 
the radiation- and matter-dominated eras and show that 
a PBH can multiply its mass by up to two orders of mag- 
nitude through the accretion of a dark matter halo. 

If PBHs can grow sufficiently by accretion (or 
if they are very large at birth) they may ac- 
count for intermediate-mass black holes (IMBHs) in 
the mass range 10 2 M Q < M < 1O 4 M , which 
have been suggested as the engines behind ultra- 
luminous x-ray sources (ULXs) recently discovered in 

1 Princeton University Department of Astrophysical Sciences, 
Peyton Hall - Ivy Lane, Princeton, NJ 08544 

Electronic address: mack@astro.princeton.edu 

2 Department of Astronomy, University of Maryland College 
Park, MD 20742 



nearby galactic disks JPewanean. Titarchuk fe Gr iffiths 



[20051 iMadhusudhan et al.ll2006UMillerll2004UMushotzkv 

2004). Based on the observed ULX luminosities of ~ 
10 39 erg s" 1 , stellar mass black holes are ruled out unless 
the emission is highly beamed. IMBHs have the appro- 
priate Eddington luminosity to explain ULXs, but since 
there is currently no easy way to produce black holes 
of this mass from stellar collap se at the abundances ob- 
served ((Frver fe Kal ogera 2001), their origins are highly 
debated. We suggest that PBHs may be able to grow 
to sufficient masses through the accretion of dark matter 
halos to account for a pervasive population of IMBHs. 

These IMBHs may also be important in the build- 
up of SMB Hs, as suggested by recent numeri- 
cal studies l|Micic. Abel fc Sigurdssonl I2005F) . a sce- 
nario that may be testable with the Laser Inter- 
ferometer Space Antenna (LISA) in coming years 
l|Micic 1 Holley-Bock elmann fc Sigurdssonl 2006| iCliffordl 
l2004HFregeau et al.ll2006() . 

Aside from observational indications, there have been 
physics-based inquiries indicating the plausibility of the 
production of PBHs through a variety of mechanisms in 
the early universe ( Carr 2005) , as briefly reviewed in |21 

One may wish to ask at this point what constraints on 
PBH production exist given current observational lim- 
its. This important question is addressed in § 18.11 and 
Figure O We will show that there is ample room for a 
population of PBHs that is both permitted and interest- 
ing. 

Our study of PBH accretion is organized as follows. In 
|2|we discuss theories of PBH formation. In JJJwe out- 
line our accretion model. Sections 01 and describe the 
accretion calculations in the radiation and matter eras, 
respectively, while in SjlD we do a combined calculation 
for both eras. In {Q we give our results for the total 
accretion possible for a PBH, and in £JEI we discuss the 
implications of our findings for the possible importance 
of PBHs in the present-day universe. Appendices lAl and 
El present details of our accretion calculations. 
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2. PRIMORDIAL BLACK HOLE PRODUCTION 

There has been a great deal of interest in primordial 
black hole production in the early universe, resulting in 
the proposal of a variety of formation mechanisms. We 
refer the reader to a review of PBHs llCarrll2005D for an 
overview of the possibilities , briefly summar ized here. In 
one mechanism of interest ( Jcdamziki ll997j) . PBHs form 
at a QCD phase transition at ~ 1M , a scale of interest 
for microlensing studies. PBHs formed at higher mass 
through other mechanisms may be natural candidates to 
solve other problems, such as that of the nature of ULX 
engines. In general, most mecha nisms creat e PBHs at 
about the horizon mass, given by llCarrll2005fl : 



general the evaporation time is given by (|Car3[2003) 
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where g e Q is the number of effective relativistic degrees 
of freedom. 

Briefly, some common mechanisms for creation of 
PBHs in the early universe are: (1) PBHs formed at 
the QCD phase transition (mentioned above), when con- 
ditions temporarily allow regions of modest overdensity 
to c ollapse into blac k holes when they enter the hori- 
zon fjedamzik I1997|) . The PBHs formed in this way 
would have a mass spectrum strongly peaked at the 
QCD epoch horizon mass (~ IMq). (2) The collapse 
of rare peaks in the density field of the early universe. 
In this case, the probability of PBH formation at a 
given epoch is determined by the nature and evolution 
of the perturbations. (3) The collapse of cosmic string 
loops ll Cgkdwgllfc Cggp er^ffiffi [G^mgaj^^gkgllar ig 
1 19931 H Hawkin d 119891: iPolnarev fc Ze mboriczl 119? 
iMacGibbon. Brandenberger fe Wichoskil 1 1998(1 . Due to 
frequent collisions and reconnections, cosmic strings may 
occasionally form loops compact enough that the loop 
is within its Schwarzschild radius in every dimension . 
(4) A soft equation of state l|Khlopov fc Pomarevl fl980'). 
If the equation of state becomes soft (e.g., during a 
phase transition), PBHs may form at peaks in den- 
sity as pressure support weakens. (5) Bubble collisions 
llCrawford fc Schrammlll982l: lHawking. Moss fc Steward 
1982: lLa fc Steinhardtl (l989 L During spontaneous sym- 
metry breaking, bubbles of broken symmetry may col- 
lide, in some cases focusing energy at a point and pro- 
ducing a black hole. In this mechanism, the PBHs would 
form at the horizon mass of the phase transition. (6) 
Colla p se of domain walls { jerezin. Kuzmin fc Tkachevl 
119831 llpser fc Sikiyiel 11984^ Closed domain walls 
forming at a second-order phase transition may col- 
lapse to form PBHs. In the case of thermal equilib- 
rium, this would result in very sm all masses, but see 
( Rubi n. Khlonov fc Sakharovl 120001) for a discussion of 
how non-equilibrium conditions may result in significant 
PBH masses. 

For cases in which the mass of the PBH is low at cre- 
ation, the PBH may evapo rate before the present day 
through Hawking radiation (:H awkinelll975|) . The limit- 
ing mass for evaporation by the present day is 10 15 g; in 
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If PBHs are to arise directly from primordial density 
perturbations, it is required that the scale of fluctuations 
set down by inflation be "blue" - i.e., the spectrum must 
have more power on small scales. In terms of inflationary 
parameters, this implies that the scalar spectral index 
n > 1, which is dis favored in the latest WMAP results 
(ISnergel et al.ll200l . 

Some mechanisms, such as the collapse of density 
peaks, may result in PBHs forming in clusters. For a dis- 
cussio n of the consequences of clustering, see ijChisholml 
I2005D. Formation via domain wall collaps e, as discussed 
in lDokuchaev. Eroshenko fc Rubml |2004) , may also lead 
to clustering, without relying on initial dark matter per- 
turbations. In that scenario, primordial black holes can 
grow through mergers to form galaxies without the help 
of initial perturbations in the dark matter. 

In this work, we assume the PBHs are rare and iso- 
lated rather than appearing in clusters, but we expect 
that clustered PBHs would increase accretion power, so 
in that respect our treatment is a conservative one. We 
may refer to a specific PBH seed mass when convenient 
for illustrative purposes, but it should be noted that our 
results are independent of the PBH formation mccha- 



3. ACCRETION MODEL 
3.1. Setup 

We model accretion of matter onto primordial black 
holes in both the radiation era and the matter era. 
In both cases we follow the calculations for ra- 
dial infall, following prev ious w ork on the g rowth 
of clusters llGunn fc Gottl Il97l lBertschingerril985l 
Fillm ore fc Goldreicbl 1198411 . Acknowledging that in a 
realistic accretion model the infall is unlikely to be per- 
fectly radial, we make the simplifying assumption that 
the angular momentum of the infalling matter causes it 
to accrete in a halo around the PBH rather than incor- 
porating itself into the PBH itself. This assumption is 
conservative from the standpoint of an estimation of the 
PBH's mass increase. 

A PBH clothed by a dark matter halo will have the 
accreting power of an object having the total mass of 
the PBH plus the halo, to the extent that the accretion 
radius (e.g., Bondi radius) is larger than the radius of 
the PBH dark matter halo. However, constraints on the 
PBH's effect on the power spectrum apply only to the 
seed mass, not to the total mass of the clothed PBH, 
since the mass accreted by the PBH is drawn from the 
surrounding matter, and the additional mass is therefore 
"compensated." In other words, a region may be defined 
around the PBH for which the overdensity is due only to 
the original PBH, with no contribution from the accreted 
mass. 

3.2. Assumptions 

In all cases, we use the cosmological para meters derived 
from the third-year WMAP data release l)Soergel et alJ 
2006). Specifically, we use the parameter set derived 
from the assumption of a flat, ACDM universe, with the 
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combination of WMAP III and all other data sets (JIa = 
0.738, n m = 0.262, h = 0.708, cr 8 = 0.751). In both the 
radiation and matter era calculations, we consider the 
accretion of dark matter only. In the radiation era, we 
assume the radiation is too stiff to accrete appreciably, as 
suggested in many pa st an alyses (iCarr fc Hawkinelll974t 
iCustodio fc Horvathlll998t iNiemever fc Jedamzikl ll9991. 
In the matter era, we ignore the small contribution to 
the PBH mass due to the accretion of baryons. 

We make the further assumptions that each PBH is 
stationary and isolated, and that the surrounding matter 
is initially in the Hubble flow. 

In all our accretion models, we end the calculation at 
a sufficiently high redshift that the effect of the cosmo- 
logical constant is negligible. For completeness, however, 
we include in Appendix [B] the outline of the calculation 
with the cosmological constant included. 

For a more detailed analysis of the consequences of gas 
accreti on onto PBHs, we refer the read er to a companion 
paper foimtti. Ostriker fc Marid lSffl^l. 

3.3. PBH velocities 

Our accretion estimate would decrease if the PBHs 
were moving quickly relative to the dark matter sur- 
rounding them; here, we assume that the PBHs are ini- 
tially stationary. We justify this assumption by consid- 
ering the likely effect of nearby density perturbations in 
the dark matter. At any epoch, we can estimate the mass 
scale at which structures are becoming nonlinear by cal- 
culating the variance of density perturbations from an 
estimate of the matter power spectrum. In Figure ^ we 
plot the nonlinear mass scales for 1 and 2a perturbations. 
The fitting formula 

M 2a = (1 x 10 17 M Q ) exp(-5.57(l + z) ' 57 ) (4) 

approximates the 2a mass perturbations. From this, 
we may calculate the characteristic circular velocity and 
thus the typical proper velocity of PBHs as a function of 
redshift: 

1/3 fl + z\ 1/2 
<V 'V : : '< 1 < km s 'i| | I —^r- I . (5) 



For redshifts down to z ~ 30, the peculiar velocities are 
low and we can consider the PBHs to be stationary. 

4. RADIATION ERA 

The details of an analytical estimate of accretion in 
the radiation era can be found in Appendix Here we 
outline the basic idea and quote the result of a numerical 
calculation. 

We begin the calculation at a redshift z w 10 7 . In 
the radiation era, the motion of a dark matter shell a 
distance r from the black hole is governed by the differ- 
ential equation 
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dt 2 r 2 

where ttibh is the black hole mass and t is time. With 
the initial conditions 



dr _ _ lr\ 
dt ~ lTl ~ 2 ti 



(7) 



at t = U, we evolve these equations forward in time until 
matter-radiation equality at z eq w 3 x 10 3 . When a shell 




Fig. 1. — Mass of 1-ct fluctuations (dotted line) and 2-<r fluctua- 
tions (dashed line) as a function of redshift. 



turns around (r = 0), we assign the matter in that shell 
to the PBH's dark matter halo. 

We find that the PBH can accrete a dark matter halo 
on order its original mass: 
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5. MATTER ERA 



The evolution of a spherically symmetric overdensity in 
the matter-dominated era has be en treated in the case of 
the growth of clusters ( Gunn & 
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glect the effect of the cosmological constant, assuming 
SI A = and fl m = 1. In the general case where A ^ 
the equation of motion of a shell of dark mat ter a radius 
r from the PBH becomes ijLahav et alJ ll9911 
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The cosmological constant term affects the accretion at 
redshifts on order 1, but since wc halt our accretions 
at higher redshift, w e fi nd growth consisten t with the 
iGunn fc Gottl <|1972|) and iBertschingerl l)1985l) result: 
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with the turnaround radius of the dark matter halo 
(which we will identify as the effective radius of the halo) 
growing with time as 



t 8 ' 9 . 



(11) 



The details of the A ^ calculation are discussed in 
Appendix IbI 

For a PBH that begins growing at matter-radiation 
equality and stops at z = 30, we find that the halo in- 
creases its total mass as 
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In the general case of a PBH growing in the matter era, 

Zeq tO Zf iS 
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6. GENERAL CASE 



In addition to approximate calculations specific to the 
matter and radiation eras, we also present the general 
result, which spans both eras and includes (for complete- 
ness) consideration of the cosmological constant. 

We start with the radial infall equation for a shell of 
matter: 



dt 2 
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For computational convenience, we recast this equation 
in terms of derivatives with respect to redshift, and we 
switch to comoving coordinates. After some algebra, we 
are left with two differential equations, one for the co- 
moving radial coordinate x(z) and one for the peculiar 
velocity of a shell v(x, z) defined by 
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v = — ktr 

dt 
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where a = l/(l + z). The integration equations take on 
a simple form in the new coordinates: 

da; — v 

d7~ir 

dv G(M acc (x) + m BH ) 

— =av-\ 

dz 
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where M acc (x) is defined as the excess matter over the 
background-density matter within the comoving radius x 
(i.e., the matter previously accreted into the halo region 
around the black hole). 

7. RESULTS 

7.1. Mass Accretion 

Our results from the combined calculation are consis- 
tent with those we obtained treating the matter and ra- 
diation eras separately. A PBH can grow by two orders 
of magnitude through the accumulation of a dark mat- 
ter halo from early in the radiation era to z ~ 30, with 
the halo mass increasing proportional to the cosmic scale 
parameter a = 1/(1 + z): 



m h (z) = 4>i 



1000 

l + z 



niBH, 



where the proportionality constant <j), 
summarizes our mass accretion result. 



(19) 

3. Figure H 



7.2. Halo Profile 

In a previous study. iBertschineerl l|1985|) performed an- 
alytical calculations of radial infall onto a central over- 
density and onto a black hole; the difference in the two 
calculations was that in the former case, the particles 
could oscillate through the center, whereas in the latter 
case they were absorbed by the black hole (as is the case 
in our simulation). Bertschinger obtained a p(r) ~ r~ 




Fig. 2. — Accreted halo mass vs. redshift. The halo radius is 
defined at an overdensity 8 = 2. We include a line to indicate the 
redshift of matter-radiation equality. 




-2.25 



Fig. 3. — Dark matter halo profile. Top panel: halo overdensity 
vs. comoving radius from PBH; bottom panel: halo mass vs. co- 
moving radius from PBH. In the inner parts of the halo, the density 
falls off as r -3 , and the profile flattens in the outer regions. 



profile for the extended overdensity and a p(r) ~ r~ 1,5 
profile for the black hole case. Our simulation resulted in 
a profile of p(r) ~ r~ 3 , differing from either of the above 
cases. Our profile is illustrated in Figure 

Since a powerlaw profile has no sharp cut-off in radius, 
we must choose a criterion by which to define the matter 
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within the halo. We may choose either the turnaround 
radius (the distance out to which shells have broken free 
of the Hubble flow) or a cut on the overdensity vs. radius; 
both criteria give similar results. Using the turnaround 
radius definition, the comoving radius of the dark matter 
halo at redshift Zf from accretion beginning at matter- 
radiation equality is 



x ta = 1.30 
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(20) 



for z f < 100. 

The radius defined by a cut on overdensity can be read 
off Figure 

7.3. Density Parameter in PBHs 

As the masses of the clothed PBHs increase, so does 
their overall density parameter. Given an initial matter 
fraction 
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the final matter fraction increases in proportion to the 
clothed PBH mass: 



UBHJ _ m BH,f 
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m BH', 



(22) 



where itlbhj includes the PBH and the accreted halo. 

In Figure |3J we illustrate that the proportional mass 
increase, while not dependent on the mass of the PBH, 
does depend on the proportion of the dark matter made 
up of PBHs. When the PBHs begin to dominate the 
dark matter, they grow less because of the decrease in 
the density of dark matter. 
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Fig. 4. — Final mass per PBH plotted as a function 
of initial mass and initial PBH density parameter, with the 
lAfshord i, McD onald &: Spergell 120031) limit included for reference 



(see : 



8. DISCUSSION & CONCLUSIONS 

Our results suggest that PBHs could grow significantly 
after their formation by acquiring a dark matter halo 



and the resulting clothed black holes could make up an 
interesting fraction of the dark matter. In the following 
discussion, we show that current observations are not in 
conflict with this conclusion, and in fact there is ample 
room both observationally and theoretically for PBHs to 
play a role in the universe today. 

8.1. Observational Limits on PBHs 

We include in Figure [5] a plot of the current observa- 
tional limits on PBHs over a wide range of masses and 
dark matter fractions. Here, we describe the limits illus- 
trated in the plot. 

Most limits on PBHs in the present-day universe are 
derived from considerations of PBHs as dark matter 
candidates. PBHs massive enough to escape evapora- 
tion would certainly qualify as "cold" and "dark" mat- 
ter; however, their existence could have noticeable ef- 
fects on processes from nucleosynthesis to galaxy for- 
mation. If PBHs form at very early times, and thus 
with very low masses (< 10 10 g), they could inter- 
fere with nucleos ynthesis by emitting part icles during 
their evaporation l)Kohri fc Yokovamal2000|) . The abun- 
dance of PBHs forming after nucleosynthesis is con- 
strained by measurements of the baryon fraction of the 
universe i|Novikov et al.llT979t) . Massive PBHs can also 
be constrained by dynamical considerations in the low- 
rcdshift universe. Some of the strongest current con- 
straints are derived from the wealth of data from mi- 
crolensing searches in Galaxy ijAlcock et al. 1 1200(1 120011 
lAfonso et al.l l2003) . The frequency and character of ob- 
served microlensing events constrains black holes in the 
mass range 0.1-1 Mq to make up less than ~ 20% of 
the dark ma tter in the Galactic halo ijAlcock et al.ll206"(i 
lOouldll2005D . The limits from microlensing are shown on 
Figure labeled "MACHO" and "EROS." 

For larger masses, constraints on PBH dark matter 
can be found by examining the effect of PBHs on the 
matter power spectrum. Calculating the excess in the 
power spec trum that Poisson-distributed PBHs w ould 
contribute, lAfshordi. McDonald k Spergell (|2003[) find 
an upper limit on present-day PBH mass of a few times 
1O 4 M , which is consistent with but tighter than the 
previous constrain t at 10 6 Mp, based on the heating of 
the Galactic disk ijLacev &: Ostrike'rlll985|h However, we 
note that this limit assumes all the dark matter is in 
PBHs: flsH = Qdm- If this assumption is relaxed, we 
find that a wide range of masses can be accommodated at 
lower density parameters. Specifically, we find that the 
product of the initial density parameter and the initial 
mass are constrained by 



m B H^BH < 



mAMS 



x 



(23) 



where x is the factor by which the PBHs increase in 
mass via accretion and tuams is ~ 10 4 AfQ. We include 
this limit in the constraint plot (Figure [HJ with the label 
"AMS." 

Other limits can be placed by considering the effects 
of compact obje cts along the l ine of si ght lensing more 
distan t sources ijWambsganssI 1200*2^ . iDalcanton et alJ 
search for the slight amplification in the contin- 



1199- 

uum emission of quasars that would be expected were 
black holes to cross the line of sight during an observa- 
tion. With a large sample of observations, they are able 



6 



to place limits on black holes in the range ~ l(r 3 Af Q to 
~ 300M©. This limit is included in the constrain ts plot 
(Figure C| and labeled "QSO." IWilkinson et~aTl ll200ll) 
place limits on the abundance of massive black holes in 
the universe based on their predicted effect of creating 
multiple images of compact radio sources. Studying a 
sample of 300 sources, they find a null result and from 
that can place a constraint on the density of black holes 
along the line of sight. Their constraint is included in 
Figure El with the label "RADIO." Coincident with the 
compact radio source study, another group found a sim- 
ilar constraint by searching fo r the same lensing ef fect in 
gamma ray burst light curves ijNemiroff et al.l200lH . The 
results of the two studies are consistent with each other, 
so for simplicity we include only the IWilkinson et alJ 
( 2001) result in Figure El 

Finally, a limit on black holes with masses of ~ WMq 
and up can be placed by observing widely orbiting binar y 
systems in the Galaxy (JYoo. Chaname fc Gouldl [2004). 
If a compact object passes between the two companion 
stars in a binary, the orbits of th e stars will be perturbed. 
lYoo. Chaname fc Gouldl ll200l use this to estimate how 
many compact objects with the ability to disturb a bi- 
nary system exist in the halo. This constraint is included 
in Figure El and labeled "WB." 

We point out that none of the above observations sig- 
nificantly constrain black holes making up less than ~ 
10% of the dark matter for a wide range of masses. 
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black hole mass (M Q ) 

Fig. 5. — Obs ervational constraints on black holes in the Galactic 
halo (see i|8.1l from microlensing experiments (EROS and MA- 
CHO), quasar variability studies (QSO), compact radio source 
lensing (RADIO), the stability of wide binaries (WB), the high- 
wavenumber matter power spectrum (AMS), and the heating of 
the Galactic disk (DISK). The region labeled "PopIII remnants" 
represents a rough estimate of the region of parameter space rel- 
evant to a scenario in which Population III star remnant IMBHs 
are responsible for ULX observations. The star symbol gives the 
posit ion of the example of a PBH ULX used in the abstract (see 
0- 



8.2. PBHs as ULXs 

This part of parameter space is consistent with an 
interpretation of the recent ULX observations as ac- 
creting, intermediate-mass black holes in nearby galax- 
ies. In this scenario, ULXs occur when IMBHs resid- 
ing in the g alaxy's halo pass through molecular clouds 
in the disk ijMapelli. Ferrara fc ReJl2006t iMillerl 12004 



I Winter. Mushotzkv fe Reynolds! 1200^) . The enhanced 
density in a molecular cloud is sufficient to trigger gas 
accretion, which causes the IMBH to emit X-ray radi- 
ation. These sources would be transient, and in any 
given galaxy the number of sources detected at any given 
time would depend on the number and distribution of the 
IMBH population and the fraction of the disk made up 
of molecular clouds. 

Several recent papers e xplore the poss i bilitie s for mak- 
ing ULXs with IMBHs. iMii fc Totanil <|2005j) estimate 
the number of ULXs expected if they are the result of 
IMBHs passing through molecular clouds. Using the 
mass and dark matter fraction estimated for IMBHs 
if the y are the compact remnants of Population III 
stars ilMadau fc Silkl 120051) . M IMBH ~ 10 2 - 10 3 and 
Qimbh /fife ~ 0.1, they find that the estimated num- 
ber of ULXs is consistent with observations. These re- 
sults do not depend on the nature of the IMBHs - if the 
IMBHs were primordial in nature rather than Pop III 
star remnants, they would also be capable of producing 
the observed sources. 

Although these results are encouraging, the present 
lack of understanding of the nature of ULXs and the 
many uncertainties that go into predictions of the con- 
sequences of a halo population of IMBHs mean that 
the issue is far from resolved. A recent paper draw- 
ing on the Mii & Totani result uses an ensemble of 
N-body simulations to place constraints on the num- 
ber of IMBHs in the Milky Way halo by drawing on 
the fact that we have not observed any ULXs in our 
galaxy l)Mapelli. Ferrara fc Real l2006f) . These authors 
find in one simulation that for a halo population of ~ 10 5 
IMBHs incorporating a fractio n 0.1% of the baryons and 
distri buted in an NFW profile l)Navarro. Frenk fc White! 
1997), the predicted number of ULXs per galaxy is on 
order 1; however, the number of lower-luminosity X-ray 
sources is overproduced. The constraint found by this 
method may be applicable to our own galaxy, but it has 
not as yet been extended to other galaxies, where ULXs 
are observed. Since the predictions depend strongly on 
not only the distribution and number of IMBHs, but also 
on the properties of the gas in the galactic disk and the 
efficiency of the black hole accretion (which in turn de- 
pends on whether or not an accretion disk is formed), it 
is difficult to generalize them to other systems. 

Other authors have suggested that ULXs may be due 
to IMBHs accreting from c aptured stellar companions 
rather than molecula r clouds llPoolev fc Rapp aportl2005l: 
iPatruno et all 120061: iMadhusudhan et al.ll2006j) . In this 
case, the ULXs would "turn on" when residing in dense 
star clusters. 

More detailed observations and simulations are re- 
quired to answer the ULX question. Here, we merely 
point out that the case for IMBH ULXs is an inter- 
esting one, easily consistent with current constraints 
on the dark matter fraction in black holes, and as we 
show, PBHs may account for or grow into IMBHs by 
the present era. Thus PBHs should be considered viable 
candidates to explain these mysterious sources. 

In Figure El we include as a region of interest the 
area in parameter s pace explored in the Mii & Totani 
l|Mii fc Totanill2005H paper (labeled "PopIII remnants" ) . 
We also mark with an asterisk the position of the scenario 
discussed in the abstract: a population of 1O 2 5 M PBHs 
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making up 0.1% of the dark matter and growing through 
accretion to incorporate 10% of the dark matter by the 
present day (we mark the mass of the seed PBH only, go- 
ing on the conservative assumption that the PBH's dark 
matter halo is not directly accreted). Both these points 
lie in a region of acceptable parameter space for the ex- 
planation of ULXs with IMBHs, but the exact extent of 
the region is difficult to define given the uncertainties 
mentioned above. 

8.3. PBHs as SMBH Seeds 

The question of the origin of supermassive black 
holes at high redshifts has attracted a great deal of 
attention since quasars have been discovered at red- 
shifts z > 6, i mplying that black holes as mas- 
sive as ~ 10 9 M^ ilFan et all l2003t iBarth et all l2003 
iWillott. McLure fc ,Iarvisll2003ft exist when the universe 
is less than 1 Gyr old. It has proven difficult to find 
a mechanism that can create such massive black holes 
so quickly. Most proposals require smaller black holes 
to act as seeds for the build-up of SMBHs. In some 
cas es, these seeds form directly from the colla pse of ha- 
los l|Lodato fc Nataraianll200alSpaans fe SirkH2006ft . In 
others, the seeds are the remnants of P opulation III stars 
l|Shapiroll2005HVolonteri fe Reesll2005j) . or they might be 
primordial. The question of whether or not SMBHs can 
be grown from the merging of stellar mass black holes has 
also been discussed in recent work (see iLi et alJ 1200(1 
and references therein). While in each scenario a case 
may be made for the ability of these seeds to result in 
the SMBHs we observe as quasars, sometimes requiring 
the i n vocation of self -interacting dark matter jQstriker 
120001 IHu et all 12006ft o r the accretion of scalar fields 
l|Bean fc Magueiiol 120021) . there is as yet no consensus 
on the matter. We suggest that PBHs may be a viable 
SMBH seed candidate because of their ability to build 



up large dark matter halos that may assist in further 
growth through baryon accretion later on. Furthermore, 
forming earlier than Pop III stars, they have more time 
to grow in the epoch before quasars are observed. 



8.4. Conclusions 

We have shown that primordial black holes can grow 
significantly after formation through the accretion of a 
dark matter halo. In the radiation era this can lead to 
an increase of total mass on order unity, while during 
matter domination, the mass can grow by roughly two 
orders of magnitude. Although a dark matter halo may 
not significantly increase the mass of the seed black hole 
itself due to the lack of a mechanism to dissipate angular 
momentum, the accumulation of a halo cannot be ignored 
when considering the ability of a PBH to accre te gas in 
later eras (in the case of collisional dark matter ijOstrikerl 
2000), the PBH can directly accrete significant amounts 
of dark matter). We have also shown that the parameter 
space available to PBHs as components of dark matter 
components is sufficient for them to play an interesting 
role in galaxies. PBHs may be viable candidates for the 
seeds of intermediate mass black holes, possibly respon- 
sible for ultra-luminous x-ray sources, or they may play a 
role in seeding supermassive black holes currently found 
in the centers of galaxies. 
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APPENDIX 

A. RADIATION ERA ANALYTICAL APPROXIMATION 
The equation describing the dynamics of a dark matter shell a distance r from a black hole of mass tubh at time t 



d 2 r Gm 
dt 2 = ~ 



Liu 

At 2 ' 



In the radiation era, we have H = l/(2t). We define r, = r(t = ij), which leads to 

dr . . lri 

where Hi is the initial Hubble parameter. 
We now consider the unperturbed solution, where uibh = 0. This reduces the above second-order equation to 

d 2 r 1 r 



dt 2 A t 2 ' 



We take the unperturbed behavior to be a powerlaw, 

ro(t) 

Then, 



' 1 i 



d 2 r _ n 
dt 2 ~ t? 



a(a-l)t a - 2 = - 



At 2 ' 



which we can solve to get 



a 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 
(A6) 
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This gives us the time evolution of the unperturbed solution: 

r 

ST* 



ro = -tV (A7) 



We can now take r = r + Sr. Differentiating twice, we get 



r + 6r- GmBH 1 r ° 1 5r (A8) 

4(. + S)' 4(2 4(! ' ' ' 

Substituting in our solutions for ro and r'o, we obtain 

GmBH f Sr 

Sr = 1 — (A9) 

The above equation is still exact, but we can take Sr/r to be small to get the lowest order solution, 

1 Sr GmBH t 
I¥ = rf~ t 
For simplicity, we now define x = Sr/ri and r = t/t,. This gives us 



Sr + ±™ I. (AID) 



^ + 7- = r- ( A11 ) 

Defining 



we have 



f X 


GuiBHtf 


4^2 " 


r 3 r 


GmBHtf 




rf 


-(*). 


1 X 

4 r 2 


e 
r 



(A12) 



(A13) 



The solution of the homogeneous equation, i + x/(4r 2 ) = 0, is x oc t°. Finding the particular integral yields 
x = — 4er, which gives a general solution 

x = At 1 ' 2 - 4er (A14) 
We solve for A by considering that at r = 1, x = 0, which makes ^4 = 4e. Thus we have 

x(r) = 4e(r 1 / 2 - r) = — (r). (A15) 

n 

For those shells of matter bound to the central black hole, each shell will have a turnaround time (time at which the 
shell ceases to expand in the Hubble flow and begins to fall back) and a collapse time (time when the radius of the 
shell goes to zero). 

We define a shell's collapse time as the time when r = 0: 

Sv 

r = r + Sr = r H r< = 0. (A16) 

n 



Rewriting this with x and r, we have 
and with the solution for x, this gives 



r i T cM + x coiiri = 0, (A17) 



1 + 4e N 2 



Tcoll = j , (A18) 

or 

tcoll = U (j-^j ■ ( Al9 ) 

To find the amount of matter accreted by the black hole in the radiation era, we choose some initial time, and find 
the amount of matter accreted between that time and matter-radiation equality. This will likely be an overestimate, 
however, as any interactions or other effects are more likely to slow accretion by pulling matter away from the black 
hole. This calculation will find the amount of matter that had sufficient time to accrete, assuming the accretion is 
steady and undisturbed. Here we perform a rough estimate of the amount of matter accreted. 

We choose as the initial time when Zi ~ fO 7 . The final time is the time of matter-radiation equality, tf = t eq , 
corresponding to z eq w 3 x 10 3 . During the radiation epoch, t oc (1 + z)~ 2 , so we have U/tf w f CT 7 . Setting tf = t co u, 
we have 

tcM = / l + 4e mm \ = iq7 ^ (A2q) 



where e m i n reflects the fact that this is a maximal estimate of the accretion. Since we defined e = Gtubh /(rft i 2 ), we 
can now write 

rl max - gr^- (A21) 



Then, 



which gives us 



4 3 ^ GrriBH n3 , . , 

lacc,max = 7; 7Tr i,maxPm,i = o^T^l Pm,o(l + z i) i (A22) 

m acc<max _ 4*Gp mfi (1 + (A23) 



With some manipulation, this becomes 



m B H 9 \t// Vi 



(1 + Zi) 3 (A24) 



where the subscript refers to the present era. This is not exact, since we are implicitly assuming that the present era 
is completely matter-dominated, i.e., £l m .o = 1. However, this can be neglected if we instead use z 30 as the final 
time - in this case, the fraction U/to is not significantly changed and the factor f2 m; o becomes unity. For e m j„, we use 
Sp/p w 2.5 x 1CT 5 . Here t f /t = 6 x 1CT 6 , t t /t f = 9 x 1(T 8 , and z t = 10 7 , which gives us 



macaco* ^ ^ (A25) 

mBH 

Thus we see that the fractional increase in mass of the black hole in the radiation era is on order 1. 

B. MATTER ERA CALCULATION WITH COSMOLOGICAL CONSTANT 

To calculate radial infall of dark matter with the effect of the cosmological constant included, we solve llLahav et alJ 
Il99l : 



d 2 r GM, Ar 



dt 2 r 2 3 ' 

The mass internal to a dark matter shell initially at ri is given by 



(Bl) 



4 

Mi = -irp„ l:i rf + m BH , (B2) 



A - = n /on-'-3 L ( B3 ) 



where Ti is the initial physical radius, with the initial density contrast defined as 

M, 
(4/3)7rf A 

From this, lLahav et alJ l|1991f) derive the equation of motion: 

^ = -~(A ci + l)A- 2 + \ i A, (B4) 

where Xi = t = Hit, and where A is the scale factor of the shell, R(t) = A(t, r,)r,. The initial conditions for the 
Hubble flow are Ai = 1, dA/dr = 1. This must be solved numerically for the initial density contrast A C i corresponding 
to the shell collapsing at the final time r/. The mass accreted from n to r/ is given by 

/ 1 — ^A,i — &BH,i \ /-n-s 

m-acc = m BH -r -p: — ~ • (B5) 



As a result of the presence of the cosmological constant, there will be a last bound shell, which is the last shell 
of matter that can in principle be accreted by the black hole. All shells internal to this are bound and will turn 
around and fall back, but tho s e bey ond it will continue expanding in the Hubble flow. Based on the calculation in 
iSubramanian. Cen fc Ostrikerl l|200ff) . we find the last bound shell to be at an initial radius of 

1 / o \ -V 3 



Z7T Pm,i \ *t>m,i, 



where ri TOj i is the initial density parameter in matter. 
The mass within the last bound shell is then 



™ A = (4Tr/3) Pmji rl, (B7) 
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and 

2 / 1 — \ 

mx = 3{-mi i ) mBH - m 

For the case of accretion beginning at matter-radiation equality, the ratio of the mass within the last bound shell to 
the initial mass of the PBH will be mx/raBH ~ 1500. This is the mass increase possible, in principle, for each PBH. 
However, as we derive this assuming each PBH to be isolated, it may occur that the PBH will run out of matter to 
accrete before this limit is reached, as it has all been accreted by neighboring black holes or that infall of much of this 
mass would occur after the present epoch. More realistically, as the PBHs grow, they will begin to interact with one 
another, and their cluster dynamics and mergers will have to be considered. The last bound shell mass can therefore 
be considered a strict upper limit on the accretion. 
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